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Abstract 

A de Sitter black hole or a black hole spacetime endowed with a positive cosmological constant 
has two Killing horizons - a black hole and a cosmological event horizon surrounding it. It is 
natural to expect that the total Bekenstein-Hawking entropy of such spacetimes should be the 
sum of the two horizons’ areas. In this work we apply the recently developed formalism using 
the Gibbons-Hawking-York boundary term and the near horizon symmetries to derive the total 
entropy of such two horizon spacetimes. We construct a suitable general geometric set up for 
general stationary axisymmetric spacetimes with two or more than two commuting Killing vector 
fields in an arbitrary spacetime dimensions. This framework helps us to deal with both the horizons 
in an equal footing. We show that in order to obtain the total entropy of such spacetimes, the near 
horizon mode functions for the diffeomorphism generating vector fields have to be restricted in a 
certain manner, compared to the single horizon spacetimes. We next discuss specific known exact 
solutions belonging to the Kerr-Newman- or the Plebanski-Demianski-de Sitter families to show 
that they fall into the category of our general framework. We end with a sketch of further possible 
extensions of this work. 

PACS: 04.70.Bw, 04.20.Jb 

Keywords: Black holes, de Sitter, cosmological event horizon, Kerr-Newman-de Sitter, Plebanski- 
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1 Introduction and motivation 

Since the discovery of the accelerated expansion of our universe, there has been intensified interest 
in the study of spacetime physics with a positive cosmological constant, A. A tiny positive A 
can very satisfactorily explain the data corresponding to the current cosmological evolution of our 
universe. Thus, black hole spacetimes endowed with a positive cosmological constant are expected 
to provide us reasonable and physically well motivated models to study the global properties of the 
black holes living in our current universe. 
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One of the most exotic features of spacetimes endowed with a positive cosmological constant 
perhaps is the existence of a cosmological event horizon, when the parameters of that solution obey 
certain conditions [I] . It can be regarded as a complementary part of the black hole event horizon, 
arising due to the repulsive effect due to positive A, at large length scales. The cosmological event 
horizon acts as a causal boundary surrounding us, as no communication along a future directed 
path is possible beyond it [1] . When we have a black hole located inside the cosmological horizon, 
we call the entire spacetime a de Sitter black hole. Thus for such black holes the natural region of 
interest is the region between the black hole and the cosmological event horizon. 

Just like the black hole event horizon, the cosmological event horizon also creates particles [I], 
and possesses thermodynamical properties, see. e.g. and references therein for some recent 

developments in this direction. 

Given that a black hole has thermodynamical properties and it creates particles, it has been 
an exciting topic to understand this from the symmetry and the microscopic point of view. In 
particular, for black hole spacetimes, the Bekenstein-Hawking entropy can solely be derived utiliz¬ 
ing the near horizon conformal properties |12j-|I8). Precisely, these formalisms use some suitable 
fall-off conditions near a horizon and investigates under action of which vector helds these near 
horizon structures are ‘preserved’. The algebra of charges corresponding to these symmetry gen¬ 
erating vector fields gives a Virasoro algebra with a central extension. Then one uses the Cardy 
formula [H]-[5D] to determine the entropy of the spacetime, which, under suitable choice of the 
mode functions, coincides with the Bekenstein-Hawking entropy. In other words, the entropy of 
a black hole spacetime is solely determined by the local symmetry at its boundary, i.e the event 
horizon. 

We further refer our reader to [21] and references therein for a recent review on various ap¬ 
proaches to understand black hole thermodynamics, the derivation of the Bekenstein-Hawking 
entropy and its quantum corrections. 

A recent approach for deriving the thermal properties of a Killing horizon can be found in [22] - 
[25] . This novel formalism uses the the Noether current associated with the variation of the Gibbons- 
Hawking-York boundary term to obtain the conserved charges associated with the diffeomorphism 
generating vector fields. It turns out that the variation of the surface term gives a boundary integral 
of the Noether current located on the event horizon. Then the requirement of the near horizon 
symmetry gives the Bekenstein-Hawking entropy of the spacetime. This method has also recently 
been applied to time dependent cosmological black hole spacetimes [26], and to conformal scalar 
hairy black holes 1221- 

Let us return to our focus - black holes in the de Sitter spacetime. Since such a spacetime has 
two Killing horizons, it is natural to expect that its entropy will be the sum of the two horizon 
areas, e.g. |2]. Our precise goal in this work is to derive this total entropy using the formalism 
of [22]-[25]. Since such spacetimes has two natural boundaries, i.e. the two horizons, we get two 
surface integrals located at the two horizons. Then we shall find out the modes of the vector fields 
that preserve both the horizons’ structure to find a Virasoro algebra ‘effectively’ representing the 
whole spacetime, giving us the total entropy (Sec. 3). To the best our knowledge, this has not been 
done before. We refer our reader to [28] for a verification of the second law of thermodynamics 
with this definition of entropy for some specific cases. See also [29] for some further discussions 
motivating this definition of the total entropy, from the point of view of area theorem. 

Let us now emphasize a rather peculiar feature of this total entropy. We consider the Schwarzschild- 
de Sitter spacetime. 
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where the smaller root is the black hole and the larger tq is the cosmological horizon. For 
3MGVA = 1, th and rc merge to l/'/A, known as the Nariai limit. The entropy and the temper¬ 
ature of the black hole and the cosmological horizons are respectively given by {Ah/4:G, kh/^tt)^ 
and {Acl4:G, Kc/27r), where A denotes the horizon area, and kh and —kc are respectively the two 
horizon’s surface gravity (with kc > 0, and kh > kc) [ I ]- 

Now, if we consider the variation of the total entropy of the spacetime, S = {Ah + Ac)/4G, we 
obtain a Smarr formula with an effective equilibrium temperature, TeS = 2tt{ih^+k,c) ’ [3E]- 

This implies that even though the two horizons have different characteristic temperatures, there 
can be an effective thermal equilibrium state when their entropies are combined. This has been 
demonstrated earlier in [30] via semiclassical tunneling method, but to the best of our knowledge, 
any field theoretic derivation of this is yet unknown. 

Thus, it is highly motivating to understand this effective thermal equilibrium state, and quite 
naturally, a first step would be to actually derive the total entropy of such spacetimes. 

The derivation of the entropy of the cosmological event horizon in any dimension can be found 
in m, where suitable fall-off condition on an asymptotic de Sitter spacetime and the near horizon 
symmetry has been used. A discussion on the relation between the Friedman equation and the 
Cardy formula can be seen in [52] ■ We also refer our reader to [55] for a discussion on phase 
transition of de Sitter black holes using the aforementioned effective thermal equilibrium state. 

The paper is organized as follows. In the next section we outline the general geometric scheme 
in which we work in. This will help us to deal with general stationary axisymmetric spacetimes with 
arbitrary number of commuting Killing vector fields in arbitrary spacetime dimensions. Precisely, 
this will provide us a timelike (non-Killing) vector field that foliates the spacetime between the 
two horizons. This vector field becomes null and Killing on both the horizons. This will enable us 
to identify Rindler geometries in the vicinity of both the horizons, and hence to treat them in an 
equal footing. Apart from the existence of the cosmological event horizon, we shall not assume any 
further explicit form or fall-off for the metric there. In other words, our method works well for the 
Nariai class de Sitter black holes as well, where the two horizons have comparable length scales. 

In Sec. 3, we use this general framework to derive the total entropy of such spacetimes extending 
the formalism of [22]-[25]. We point out there that the mode functions corresponding to the 
symmetry generating vector fields near the horizons have to be restricted in a certain manner, 
compared to the single horizon spacetimes, in order to obtain the total entropy. Sec. 4 is devoted 
to address known non-trivial explicit examples from four and higher dimensions, to demonstrate 
they all fall under the general framework we built. We also address two cases of non-minimal 
couplings here and show in particular, some qualitative difference for the Brans-Dicke field for our 
case, when compared to the asymptotically flat spacetimes. Finally, we discuss our results in Sec. 5. 

We shall work with mostly positive signature of the metric and set c = fcs = = 1 throughout, 

but will retain Newton’s constant, G. In different spacetime dimensions, different values of it will 
be understood. 

2 The general near horizon geometry 

We shall derive below the general geometric framework we will be working in. The first part of which 
essentially deals with the construction of Killing horizons in stationary axisymmetric spacetimes 
of general dimensions with two or more commuting Killing vector fields. This will help us to deal 
with both the black hole and the cosmological horizon in an equal footing and in a much convenient 
manner than dealing with exact solutions case by case. The essential details of this can be found in 
e.g. [HU [35] and references therein. For the sake of self consistency and convenience of the reader, 
we shall briefly outline them here. 

The next part consists of identifying an (1 -|- l)-dimensional geometry in a general way, such 
that it becomes the Rindler on any of the two Killing horizons. 

We assume that the spacetime is an n-dimensional torsion-free manifold and satisfies Einstein’s 
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equations. We assume that the spacetime is stationary, axisymmetric and is endowed with two or 
more than two commuting Killing vector fields, 

V(a4) = 0 = 

= 0 = {i,j = l,2...m, with m < n) (3) 

where and respectively generates stationarity and axisymmetries and the vanishing com¬ 
mutators are represented by the vanishing Lie derivatives. Since the spacetime is stationary, not 
static, we take ^ 0 for alH = 1 , 2 ,... m. Also, to allow sufficient generality in our method, 

we further assume that the axisymmetric Killing vector fields are not mutually orthogonal 

as well. There can be additional spatial isometries orthogonal to the stationary Killing vector field 
but for our present purpose we need to worry only about isometries non-orthogonal to 
For convenience, we shall first discuss the case of three commuting Killing vector fields (^, 
Generalization to higher numbers or specialization to two such commuting fields will be clear from 
this, as we shall see below. 

We assume that the (n—3)-dimensional spacelike surfaces orthogonal to these three commuting 
Killing vector fields form integral submanifolds, which essentially means the vector fields span¬ 
ning the subspace form a Lie algebra between themselves, which in turn implies Frobenius-like 
conditions |35) . 

= 0. (4) 


Clearly, the chief difference between the static and stationary axisymmetric spacetime is that, for 
the later the timelike Killing vector field is not hypersurface orthogonal. For our convenience, we 
shall now construct a foliation of the spacetime by constructing a hypersurface orthogonal (non- 
Killing) vector field. To do this, we define a 1-form Xa as 


Xa = Ca + + a2(x)(^i^\ 

so that = 0 = identically everywhere, giving 


ai{x) = 




0 : 2 ( 2 ;) = 




(^(I).^(2))2_j2y2 > (^(1).^(2))2_J2J2 

where we have written for the norms, = +ff and ■ (jP''> = -l-/|. Let = 

Then the norm of Xo is given by 


(5) 

( 6 ) 

-A2. 


xax“ = -P=\ -y + 




(^( 1 ). 0 ( 2 )) 2 _^ 2 ^, 


(7) 


Since for any two spacelike vector fields A and B, we always have A - B < |A||i3|, the denominator 
of the second term on the right hand side is negative. This shows that Xa is timelike when /3^ > 0. 
The price we have paid doing this orthogonalization is that, Xo is not a Killing field in general, 

V (aXb) = <(>(a Vb) Oi (x) -b V 6) 02 (x) . (8) 

In terms of Xa, the first of Eq. (jd)) can be written as Xc^dXe] = 0, which admits general 

solution 

V[aX&] = ^^l[aXb] + + i^lX[a<^i[^ + i"2X[a</>if + , (9) 

where (z = 1,2,3) are 1-forms orthogonal to x“i and and Vi[x) {i = 1,2,3) are 

functions. These functions and 1-forms can be determined exactly, chiefly using the commutativity 
of the Killing vector fields [34] . 

V[aXb] = {Xb^aP^ - Xa^bP^) , (10) 
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which shows Xa satisfies the Frobenius condition, X[a'^bXc] = 0; hence is orthogonal to the 

family of (n — l)-dimensional spacelike hypersurfaces containing t/fa^’s. This is a crucial result for 
proceeding further. 

Having obtained the foliation of the spacetime, we now proceed to define the Killing horizons. 
Eq. (flUl) shows by the torsion-free condition that, for any 0^=0 hypersurface (say "H), 

X[b^a]/3 |^2_^q P I ^2_,,q ^ 0; (H) 

so that on any such hypersurface 'H, we may write 

^ 2K{x)Xa, (12) 

where k{x) is a smooth function defined on Ti.. 

The next step is to prove that any such compact surface "H is a Killing horizon, in the sense 
that the functions ai{x) (Eq.s (O, ([5])) becomes constant on T-L. This involves constructing a null 
geodesic congruence for ka = (with x°'^aT '■= 1), and then solving for the Raychaudhuri 

equation on TL. We shall not go into the details of this proof here referring our reader to [33] and 
references therein for this. 

Then following similar steps as in four spacetime dimensions [35] . we can show that k is a 
constant on 'H, and is given by 


2 (Va/32) (V“/32) 

4/32 


(13) 


known as the surface gravity of the Killing horizon. We shall assume in the following k ^ 0, always. 

Thus we have seen that the foliation timelike vector field x“, smoothly becomes the horizon 
Killing vector field (say, Xn)- Sitter black hole spacetimes we wish to deal with, we have 

two such compact /3^ = 0 surfaces. The smaller one is the black hole event horizon and the larger 
one is the cosmological event horizon, and the vector field smoothly becomes null and Killing 
on both of them. 

The next step is to show using this general framework that, we can select a part of the near 
horizon geometry, which is Rindler-like. In order to see this, we define a 1-form 

Xa := -Va/3, (14) 

n 

which is orthogonal to and as can be seen by using ©■ Also, Eq. (IT31) shows, XaX°' = -|-1. 
Let X be the parameter along X“, such that X“VaX := 1. We have from the action of a vector 
field on functions [55] . 


XaX“ = l=iA“Va^=-^, (15) 

K KdX 

which gives, X = ^. If we choose X“ to be one of the basis vectors orthogonal to x“j it is clear 
that the metric infinitesimally close to H takes the form 

9ab — ~ X'HaX'Hb T -\- Xab: (^^5) 

where the spacelike compact (n — 2)-section jab is orthogonal to both x“ and X“. It is clear that 
while xh is tangent to H, the vector field X“ defines orthogonality or ‘away from’ H. The ‘x — X’ 
part of the near horizon coincides with the Rindler metric. 

We assume that the basis vectors spanning jab have neither vanishing nor diverging norms. This 
is just because otherwise we will have either vanishing or diverging horizon ‘area’ (:= J (detqah)). 
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We note that the function k^{x) = ^ ^ , where /3^ is not necessarily vanishing this time, 

smoothly coincides with (Eq. (TT^ l on 'H. Then it is clear that with the vector field 
we may write the general spacetime metric as 

gab = —P "^XaXb + ^aXb + ^ab, (17) 

which smoothly coincides with (IT6l) on any of the two Killing horizons. This helps us to deal with 
the black hole and the cosmological horizon in an equal footing. 

To summarize, we have found a foliation of an n-dimensional stationary axisymmetric spacetime 
with three commuting Killing vector fields along a timelike vector field. Whenever that vector field 
becomes null on a compact surface, it becomes Killing as well, making the null surface a Killing 
horizon. For our concern, we have two such Killing horizons. We have also shown that the general 
spacetime metric (HZ} coincides with (flfill . in an infinitesimal neighborhood of any of the horizons. 

For two commuting Killing vector fields, we set any one of and (and hence / 12 , Eq.s ([6])) 
to zero. Also, there is only two conditions analogous to dH) now. They include ^ and any one of 
the axisymmetric Killing vector fields. 

For more than three commuting Killing vector fields, say four, we assume that the (n — 4)- 
subspace orthogonal to those Killing vector fields form integral submanifolds. Accordingly, we have 
four Frobenius-like conditions dH). We proceed then as earlier to define Xa as Xa = Ca + with 

i = 1,2,3. We find ai by imposing the orthogonality between Xa and the axisymmetric Killing 
vector fields and proceed as earlier for the rest of the construction. This process may go on and 
clearly can accommodate arbitrary number of commuting Killing vector fields. 

For our convenience we further define a new coordinate X = JkX /2 in (IT6l) to get 


9ah 2^j^X'HaX'Hb T {‘^bvX'jXfj^Xb Xabi 


(18) 


where A“ is the tangent vector field associated with the new coordinate A, and AqA“ = (2kA)“^. 

Finally, we note that since the Killing vector fields commute, we may specify coordinates along 
them, at least locally. On the other hand, since the horizon Killing vector field Xu ^ linear 
combination of those Killing fields with constant coefficients, it is clear that we can treat Xu 
a coordinate Killing vector field. We shall denote the surface gravities (k) of the black hole and 
the cosmological event horizon by kh and —kc (with kc > 0) respectively and will always work 
with the absolute value of the cosmological event horizon’s surface gravity, in order to maintain the 
correct signature of the metric (flSl) . 

With all these necessary geometric ingredients, we are now ready to go into the derivation of 
entropy. 


3 General derivation of the entropy 

We shall use below the formalism developed in [12]-[2n] using the Gibbons-Hawking-York surface- 
counterterm in order to calculate the total entropy of a stationary axisymmetric de Sitter black 
hole spacetime. 

Let us briefly review the formalism first. The Gibbons-Hawking-York surface term subject to 
the variation at the boundary is given by, 

Asu. = ^[ [d--^x]K=^ [ [d-x]XUKix)N-ix)), (19) 

JdM ottG J 

where G is the Newton constant in dimension n, and \dr'~^x\, [d^x] stand for the invariant volume 
measures in respective dimensions. K is the trace of the extrinsic curvature of an (n—l)-dimensional 
boundary hypersurface, dAi . In the second integral, which is valid on the entire spacetime manifold 


6 







Ai, K{x) and N‘^{x) are respectively a function and vector field that smoothly coincide with K 
and N°' on dAi, where N°‘ is the unit normal to dAi. 

The conserved Noether charge corresponding to the variation of the integrand of the second 
integral of (1191) under infinitesimal diffeomorphism generated by a vector field is given by 

QIC] = \j = ^J (20) 

where J“ is the conserved Noether current, is an antisymmetric tensor field given by, = 
and is interpreted as the Noether potential, as T“[C] = S is a suitable 

hypersurface, and the second integral is the (n — 2)-dimensional boundary of S. The choice of 
S is made in such a way that its boundary coincides with the (n — 2 )-section of dAi. 1/7 is the 
determinant of the induced metric on that boundary, and ^/^d'yab = — x{NaAIb — NbAIa) 

is the area element. The vector fields N°' and are chosen to be unit spacelike and timelike, 
respectively. 

For a black hole spacetime, the natural choice of the (n — l)-dimensional hypersurface in (fT^ 
is clearly the event horizon. The hypersurface E in (I20I) is a spatial hypersurface, and the {n — 2)- 
dimensional subspace in (I20I1 is the compact spatial section of the event horizon, spanned by angular 
coordinates. 

The bracket algebra of charges o generated by different vector fields is given by 

[QlCmhQlCn]] := i(^uQ[Cn] - <5c„Q[Cm]) = ^ ^ y^d^abiCAU ” C?nACn]), (21) 


where (fc,„Q[Cn] := J^dEaAcmiVdet gJ°-[Cn])- 

The next step is to identify an infinite discrete set of diffeomorphism generating vector helds 
{Cm} which leave the near horizon geometry invariant. It can then be shown that for such vector 
fields, (j21l) can be identified with the Virasoro algebra with a central extension. One then uses the 
Cardy formula [Bun] in order to compute the entropy of the spacetime. 

For the de Sitter black holes, we have two natural boundaries - the black hole horizon along 
with the cosmological event horizon. In other words, the so called ‘bulk’ of the de Sitter black hole 
spacetimes is the region between these two horizons. For such two natural boundaries, the first 
integral in Eq. (HI splits into two pieces - on the two hypersurfaces located at the two horizons. 

Accordingly, the Noether charge in (1201) will consist of two integrals at the two Killing horizons, 
similarly for the algebra satished by the charges, Eq. (EH). The hypersurface S is orthogonal to 
the foliation vector field derived in the previous section. Since smoothly becomes null and 
Killing {xu) on both the horizons, we have obtained the two boundary integrals in a quite natural 
manner. 

Let us first evaluate the charge corresponding to the horizon Killing vector field Xu ^9- (EO])- 
We have seen in the previous section that the spacetime metric formally looks the same in the 
neighborhood of both the Killing horizons, m- We choose for the black hole horizon, = 

The trace of the extrinsic curvature at the black hole event 


Xu 


•\/2kh X 


^2khXX'^, and = 

horizon is given by Kjj = — Likewise, we get the unit vectors and the trace of the 

extrinsic curvature on the cosmological horizon, by replacing kh with kq- We hnd 


Q = IJ dXar = i ^ Vidxabr^ + IJ^ Vidxabr^ = 


khAh + kcAc 

’ 


( 22 ) 


where we have defined the (n — 2)-dimensional ‘area’ as A = f ,yxd"~^x, which correspond to the 
compact spatial sections at the two horizons. 

We shall next obtain a Virasoro algebra for the charges generated by vector fields which preserves 
the near horizon structures. Eq. (EH for our case becomes. 


[Q[Cm], Q[Cn]] = o / AdlabiCACm] - C^aACn]) + T: / VldjabiCACm] - Cm^lCn])- (23) 


IH 


'C 
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We shall look for the set of vector fields {C“}, which has only non-vanishing ‘time’ and X compo¬ 
nents, with respect the spacetime metric (1X51) . Let us collectively denote the spatial coordinates 
and basis vectors tangent to the horizon by {0*} and {0“} respectively, spanning in Eq. (ITSl) . 
Let 3>(a:) be the norm of any of the basis vectors {0“}- Clearly, in order to have the area of the 
horizon to be finite and non-vanishing, $ has to finite and non-vanishing, too. Thus, in the neigh¬ 
borhood of a Killing horizon, we may expand $ = <f>i(0*) -t- 0{fi) + higher order terms in /3. Also, 
since 0“’s are tangent to the horizon, where /3 = 0, we must have 0“Va/3 —?> 0 in the infinitesimal 
neighborhood of any Killing horizon. Putting these all in together, using Eq.s (nsi), (nnD, (dll) and 
the chain rule for the partial derivatives, we get 

dx^ = {dp^){d~l3){dxX), (24) 

to be vanishing in the neighborhood of any Killing horizon. This is analogous to the static and 
spherically symmetric case : the metric functions spanning the 2-sphere (at r = rn) do not depend 
upon the spacelike Rindler coordinate. A, on or in the infinitesimal neighborhood of the horizon. 

In order to preserve the near horizon geometry, we must impose £c,gTT = 0 = £c,gxx, where 
r represents the parameter or coordinate along the horizon Killing vector field Xk- Solving these 
two equations involves only the Rindler part (r — X) of the metric, and one obtains [22l 123) . 

r =r(r,A,0*)-^a,T(r,A,0O, =-Aa,T(r, A, 0^, (25) 

where T is some smooth but otherwise arbitrary function. In terms of this vector field, Eq. (1221) 
reads 

« = Sg /„ + Sg X 

We now expand the function T in terms of infinite number of discrete eigenmodes as, T = ^ AmTm 

m 

with m integer, so that for each m, we call the corresponding vector field as Cm- It is usual to 
choose Tm = where (j/^s are the parameters along the axisymmetric Killing 

vector fields, Iq and l£s are constants, g{x) = —Iq f With this choice of the modes, Cm’s 

satisfy an infinite dimensional discrete Lie algebra over a circle, 

[Cm^ Cn]ljB ~ i(m 7l)Cm-|-ni (27) 


where the subscript ‘LB’ denotes the Lie bracket. 

Thus the modes formally look the same on both the horizons, but the eigenvalue Iq may be 
different. We shall call them as Iqh and Iqc respectively, for the black hole and the cosmological 
event horizon. Moreover, we have to fix them uniquely as well in order to derive the entropy of the 
whole system, as we shall see below. This is qualitatively different from the single horizon system 
discussed in [22l [23] , where one can leave lo completely arbitrary. 

Now, due to the axisymmetric geometry, the mode functions must be periodic in the Killing 
parameters (j£ of the axisymmetric Killing vector fields. Then, since we have assumed the horizons 
to be compact, (j£'s are tangent to them and (1^ becomes 


Q 


m 


{kh/Ioh)Ah + {Kc/ioc)^C 


^m,0- 


Likewise, the algebra of the charges, Eq. (ED) gives for our two horizon spacetimes. 


(28) 


Qn\ — 


i{m — n) 
SttG 


{{kh/Ioh)Ah + {kc/Ioc)^c) ^ m + n ,0 


T— {Ah{Ioh/kh) + Ac(Zoc/kc)) 6m+n,0, 

IdttG 


( 29 ) 







which is a Virasoro algebra, effectively encompassing both the boundaries, and hence the bulk. We 
can identify the zero mode energy or the Hamiltonian and the central charge, C from Eq.s (I28L 

(ESI), 


p _ {kh/Ioh)Ah + {kc/Ioc)Ac C Ah{Ioh/kh) + Ac{Ioc/kc) 

~ ’ 12 “ leTrG 


(30) 


According to the Cardy formula [El |20] , the entropy of the system is given by S' = 27r y . 
Then it is clear that in order to get the entropy, we must set Ioh = kh and Iqc = in (P|. This 
choice gives the same modes as one obtains via the method of ‘asymptotic’ fall-off near the horizon, 

e.g. m]. 

With this choice, we get via the Cardy formula 


S = 27r 



Ah + Ac 
4G 


(31) 


i.e. the total Bekenstein-Hawking entropy of the de Sitter black hole spacetimes. 

Before proceeding further, let us summarize what we have done so far. Since the de Sitter black 
holes have two Killing horizons, the Gibbons-Hawking-York surface counter-term splits into two 
pieces, corresponding to the two horizons. The hypersurface E appearing in Eq. (12011 is the one 
which is orthogonal to the timelike vector field y“, derived in Sec. 2. In terms of the two boundary 
integrals and choice of appropriate mode functions on the horizons, we have actually derived the 
total entropy of such spacetimes. 

We have also seen that as long as the derivation of the total entropy is concerned, the choice 
of the mode functions are much more restricted than the single horizon spacetimes. Clearly, doing 
so is absolutely justified. We start with considering vector fields generating diffeomorphism in the 
entire spacetime, and due to the existence of boundaries, we only consider their explicit forms on 
the boundaries themselves. Since our spacetime has two horizons, the diffeomorphism generating 
vector fields Cm’s assume different forms on them. In other words, this analysis can be regarded as 
‘doubly local’ instead of ‘local’ [53 [E], as the single horizon systems. Most importantly, both the 
surface gravities can be formally expressed by m, so that kh and kc are nothing but the /3^ —>■ 0 
limits of the smooth function k(x) = y/ (Va/3^)(V“/3^)/4/32. Then it is clear that even though Iqh 
and Iqc have different values numerically, they are formally exactly the same. 

We shall consider below some explicit and non-trivial exact solutions in order to demonstrate 
that they indeed fall under the scope of the general analysis we have done so far. 


4 Explicit examples 

4.1 The Kerr-Newman- and the Plebanski-Demianski-de Sitter families 


Let us begin with the Kerr-Newman-de Sitter spacetime in four spacetime dimensions, whose metric 
in the Boyer-Lindquist coordinate reads 


where 


ds'^ = — 


Ar — sin^ ffAg 


2 a sin 0 


dt^ -— ((r^ -|- a^) Ag — A^) dtdcj) 


+ —5^ ((r^ + Ag - Ara^ sin^ 0') dtf^ + ^dr'^ -|- 


d0^ 


Ar = {r^+a^){l-Kr^/^)-2MGr + q^, Ag = 1-f Aa^ cos^ 0/3 

S = 1 -f Aa^/3, =r‘^ + cos^ 0, 


(32) 


( 33 ) 
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with M, a and q are respectively the mass, rotation parameter and charge. For a Dyonic black 
hole, is understood as the sum of the square of the electric and magnetic charges. 

Firstly, since the timelike and axisymmetric Killing vector fields are coordinate helds = 
4>°' = they trivially commute^. The 2-planes orthogonal to these Killing vectors are 

spanned by the coordinate vector fields, (9^)“ and Hence they also commute to give a trivial 

Lie algebra and thus form integral submanifolds [35]. This was a crucial assumption made in Sec. 2. 
Also, as we have seen in Sec. 2, this guarantees the existence of the hypersurface orthogonal timelike 
vector field = ^ - (? • </'/«<> • = df - 

The black hole and the cosmological event horizon of (l32|) correspond to the largest and the 
next to largest positive roots of A^. = 0. The inner or the Cauchy horizon will not concern us for 
our present purpose. 

The norm —of x“, close to any of the horizons (A,. —0) is given by 



g4>4> 


(r2 -h a^f 


0{Al), 


(34) 


which is null on both the horizons. Also, it is eassy to check from the metric functions ([32]) that 
the function a = 4^ becomes a constant whenever A^ = 0. 

Thus the Kerr-Newman-de Sitter spacetime falls into the general geometric category we dis¬ 
cussed in Sec. 2. 

Let us now discuss the derivation of the Rindler coordinate following Sec. 2. According to 
Eq. (11411 . infinitesimally close to any horizon, we choose 




Vg/? 

Kff.c’ 


(35) 


where the subscripts or superscripts ‘H,C’ denote black hole an the cosmological horizon respec¬ 
tively. Using Eq.s (15^ . (IMl) . we have {dg^Xa ~ in the infinitesimal vicinity of any 

of the horizons. Likewise, for the Rindler coordinate (Eq. (IT^ l. it is easy to check that 
{de)°‘Xa ^ (!l(Ay^). Thus infinitesimally close to the horizon, the vector fields dg becomes orthog¬ 
onal to Xa or A“. 

So, the Kerr-Newman-de Sitter spacetime takes the form given in Eq. (TTSl) infinitesimally close 
to the horizons with jab spaanned by {dg)°' and (9^)“. Thus, it falls into the general scheme 
discussed in Sec.s 2 and 3. 

An exact asymptotically anti-de Sitter black hole solution with two independent rotation pa¬ 
rameters in five dimensional minimal supergravity can be seen in [36] . Erom this solution we can 
obatain a de Sitter black hole spacetime via analytic continuation |10j . 


ds^ = - 


Ag (l — Ag {2Mp^ — q^ — ‘^abqg'^p^) 


-‘a'—'b 




de + + ^do^ 

Ar Ag 


(e a?) sin^ 6 (2MGp^ — q^) sin"^ 9 + 2abqp^ sin^ 9 


{e e) cos^ 9 e (2MGp2 _ 0 _|_ 2abqp^ cos^ 9 




d(je 

de^ 


2Ag sin^ 9 [a {2MGp^ — q^) + bqp^ (l — 


P 


2Ag cos^ 9 [b {2MGp^ — q^) -|- aqp^ (l — 


P 


dtdcj) 


dtdip 


^Note that in our general analysis, we did not need to assume that the Killing fields are global coordinate fields 
everywhere inside the bulk. 
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where 


2 sin^ 6 cos^ 9 \ah [2MGp^ — q^) + qp^ (a^ + 6^)1 
+ - „ - d(t)d'il}, 


P'^^a^b 


( 36 ) 


p^ = [r^ + cos^ 6 + sin^ 9) , Ag = (^1 + cos^ 9 + sin^ 9) , = (1 

(r^ + a?){r^ + b^){l — g'^r'^) + q^ + 2abq 


— (1 + ^ ff )i — 


- 2MG 


0^9^). 

(37) 


The parameters M, a, b, q specify respectively the mass, two independent rotations and the charge 
of the black hole and g^ is the positive cosmological constant. 

Clearly, being coordinate fields, all the three Killing vector fields (9t)“i {d^Y commute. 

The 2-planes orthogonal to them are spanned by coordinate fields along r and 0, which commute 
and trivially form a Lie algebra. Thus those two planes are integral submanifolds. This gurantees 
the existence of the timelike vector field orthogonal to the (r, 0, (j), Y) family of hypersurfaces, 
which from the discussions of Sec. 2 is written as 

,,a /p,\a {9tri>giptp — gttpgtfjTp) f ^ {gtijjg^i^ — gt^g4iijj) f^ ^,a 

X =[Ot) -7 - 7 - -7 - 7 -■ ('5^) 

The horizons of this spacetime are the positive roots of = 0. A discussion on thermodynamics 
of this spacetime can be seen in m- Infinitesimaly close to any of the horizons, the norm of 
takes the form. 


X Xa = 


p^r'^Ar 


[{r^ -I- aY{r^ + 0^) + abqY 


0{Al). 


(39) 


Thus x“ becomes null on the horizons. It is easy to check from the metric functions that infinites¬ 
imally close to any of the horizons, coefficient functions in (1381) become constants 


i.94>4'9i’'4’ ~ (Sb^)^) 


r=rH, rc 


{gt-pgrp^} gt4>g4nij) 

{94>4>9ip4> ~ 


a{r'^ -I- 6^)(1 — 5 ^r^) -|- bq 
{r^ + aY{r^ + bY + abq 


r—TH ‘^C 


b{r^ -I- Y){1 — g^r^) + aq 
(r^ -|- a^)(r2 -|- bY + abq 


r=rH, rc 


(40) 


where th and rc denote respectively, the balck hole and the cosmological horizon radii. 

As earlier, it can also be checked easily that the near horizon Rindler coordinates can be defined 
and 7 a 6 in Eq. m for this case is spanned by (0, </>, Y)- Thus the general analysis of Sec.s 2 and 3 
holds perfectly for (l36|) . 

We shall next consider the Plebansky-Demianski-de Sitter black hole spacetimes. This class is 
a geneailazation of the Kerr-Newman family, in the sense that apart from mass, charge, rotation 
and the cosmological constant, it contains additional parameters. The complete family of the 
Plebanski-Demianski-de Sitter class spacetimes which might represent de Sitter black holes has the 
metric [57] . 




Ar 


2 2 
P 


P 


— ^ idt — i asin^ 0 -|- 4^ sin^ ^ j j -|- ^^dr^ + ^ {adt — -I- (a -I- /)^) -I- ^ sin^ 9d9 


P 


( 41 ) 


where 


fl = 1 — — (Z-I-acos0) r, p^ = r'^ + (I + acosdY , P = sin^ 0 (l — 03 cos0 — 04 cos^ 0) 

UJ ^ ' 

Ar = (w^fc-I--I-- 2MGr-I-er^ - -I-A/3) (42) 


ijj 
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The parameters a, w, q, q^, e and k are arbitrary, and and 04 are determined via them. These 
parameters have their physical meaning in certain special sub-classes only. 

For example, for a = 0, the above metric becomes the Kerr-Newman-NUT-de Sitter solu¬ 
tion [37] . 

ds^ = —[dt — (asin^ 0 + 4,1 sin^ 6/2)d(j)] ^ -I- "I —2 ~ + (a + lY)d(f\ ^ ^ sin^ (43) 

where 

= r'^ + {I + acosO)^ , P = sin^ 0^1-1- cos 9 + 

Ar = {a^-l^ + q^ + ql,)-2MGr + r^-A{{a^-l^)l^ + {a^/3 + 2fy+r'^/3), (44) 



where q and q^, are electric and magnetic charges and I is the NUT parameter. 

We shall consider the most general class given by gu, implicitly assuming it represents de Sitter 
black holes. The black hole and the cosmological horizons are the two largest roots of = 0. It is 
easy to argue the existence of the hypersurface orthogonal timelike vector field x“ as earlier. The 
norm of x“ behaves as A,. —>■ 0 as, 

/32 = [^2 ^ ^ 1^2 _ gjj^2 0/2)2] + 0(A2), (45) 

which is null. It can also be verified from the metric functions that the function {£,-(j ))/becomes 
a constant on A,. = 0. Thus the most general class of Plebanski-Demianski metrics (HD) falls into 
our geometrical assumptions, and clearly, when it represents a de Sitter black hole spacetime (such 
as ra), we can derive its total entropy. 

Finally, we shall briefly discuss the Kerr-de Sitter spacetime in generic spacetime dimensions 
n |38l |39] . The metric of which reads in the Boyer-Lindquist like coordinates [ini IMl [39] , 


ds^ = -W{l- g\^) d^ + ^ Iwdt-Y,'^ 


ainfdcjii 


U 


^ (r2-|-aj)2 / 2 j.2 , J 2\ Udr'^ 

2^ md4>t + dpt) 


+er^dv‘^ + 


IU(I -g2r2) 


X-2MG 
Pidpi + er'^vdiA , (46) 


N 


where 


2A 

U 


N 


N 


{n-l){n-2)g^ W = Y^gy^. + ev\ A = (l - 5 V) 


2=1 


2=1 


z 


N 


1 - 


i-E 


2,,2 




li = l + g'^a^, 


(47) 


i=l 


where N is the integer part of (n — l)/2. The constant e is -l-I (0) for even (odd) spacetime di¬ 
mensions. (/)i’s are the coordinates along the axisymmetric Killing vector fields, and Oi are the 
corresponding independent rotation parameters and g^ is the cosmological constant. The coordi¬ 
nates Pi and v are not independent, but are related via the constraint. 


N 

2=1 


= 1 , 


(48) 


ensuring the correct dimensionality of the spacetime. 
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Discussion of thermodynamic properties of this spacetime including that of the variation of the 
cosmological constant can be found in m- 

The spacetime is endowed with total (N+l) Killing vector fields. Since all of them are coordinate 
Killing vector fields, they commute. The {n — N— l)-dimensional subspace is spanned by coordinate 
vector fields (3;,)“}. Since the spatial coordinates along these vector fields satisfy the 

constraint (I48L we may replace by using the chain rule of the partial derivatives, any specific dfj,. 
as a linear combination of partial derivatives of the set d,^}, j ^ i- Then it is clear that the 

Lie brackets like or [5^., always be a linear combination of the coordinate 

vector helds spanning this subspace. This is a Lie algebra, but not trivial as the earlier examples. 
Nevertheless, this ensures that the subspace orthogonal to the Killing vector fields form an integral 
submanifolds [35]. 

Then from this, the existence of the hypersurface orthogonal vector field x“, and its Killing 
property when it becomes null follows as earlier. 


4.2 Non-minimal couplings 

Before we end, we shall address two cases of non-minimal couplings, the hairy black hole with 
conformal scalar held and the Brans-Dicke theory in the Jordan frame. Since for any non-minimal 
coupling the Ricci scalar term in the action gets modihed as, f{(p) R, where ip is the scalar held, 
the presence of such coupling modihes the surface counterterm in Eq. m, as well. 

For static and spherically symmetric de Sitter black hole spacetimes, a derivation of Bekenstein- 
Hawking-Wald entropy using the surface counterterm has recently been done in for the black 
hole event horizon. The surface term reads in this case, 


A 


sur 



27rG(n - 2)ip'^ \ 

(n - 1) ) ' 


(49) 


We shall briehy discuss this case for stationary axisymmetric spacetimes. It was shown in m that 
a black hole with conformal scalar hair and a positive cosmological constant cannot have a slow 
rotation in four spacetime dimensions. As was also argued there, it may be possible though, that a 
solution with generic rotation exists - at least one cannot rule out the possibility. In the following, 
we shall assume that such a solution indeed exists, and it falls into the geometric category described 
in Sec. 2. 

Following the conserved Noether charge Q corresponding to the diffeomorphism, instead 
of Eq. (1251) . in this case becomes 


Q 



2'KG{n — 2)ip^ \ 

i.n-1) ) 

27rG(n — 2)(/3^ \ 

(n - I) J 


(KcT-iaj). 


(50) 


We note that for axisymmetric spacetimes, the field p may depend on the non-Killing coordi¬ 
nates tangent to the horizons (such as the polar angle 0). Thus unlike the spherically symmetric 
spacetimes [27], we cannot pull out the scalar field out of the integration in Eq. (|50l) . 

Using the suitable mode decomposition as described in Sec. 3, we find 


— 


SttG 


1 - 


2TTG(n - 2){ip^)H 
(n- I) 


{khAh)/Ioh + ( 1 “ 


27rG(n - 2){p^)c 
(n- 1) 


(kcAc)/Ioc 


^771,0? 


where on any of the horizons we have defined. 




1 

Ah,c 





(52) 


(51) 
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Likewise, we find the algebra of charges 


[Qm: Qn\ — 


i{m — n) 

8ttG 

im^ 

IQttG 


1 - 


1 - 


27rG(n - 2)(y>^)g 
(n- 1) 

27rG(n - 2)((/?^)g 
(n- 1) 


{kh/Ioh)Ah + 
AniloH / kh) + 


1 - I (.a/kc)Ac 


1 - 


(n- 1) 
27rG(n - 2)(v?^)c 
(n- 1) 


Ac{Ioc/kc) 


^m+n,Q 

^m+n,0- 


( 53 ) 


Setting Iqh = kh and Iqc = kc, obtain the entropy. 

Finally, we shall discuss thermodynamics of black holes in the Brans-Dicke theory (see, e.g. [JT] 
and references therein, for discussions on asymptotically flat spacetimes). The action of the 
Einstein-Brans-Dicke theory with a cosmological constant in the Jordan frame reads (see e.g. [42] 
and references therein) 


S = J [d^x] 


(pR-2A-- 


(54) 


where p is the Brans-Dicke scalar field, and lo is called the Brans-Dicke parameter. The inverse 
of the Brans-Dicke field, acts as a spacetime dependent or dynamical Newton’s ‘constant’. 
Also, for w = oo, the field becomes a constant and the theory coincides with the Einstein gravity, 

^ = =T^- 

The derivation of the surface couterterm for the Brans-Dicke theory is similar to that of (14^ . 
giving 


= 2 [ [d^-^x]pK. (55) 

Now, for stationary asymptotically flat black hole spacetimes in four spacetime dimensions, the 
Brans-Dicke theory obeys a no hair theorem [53], which states that p is necessarily a constant in 
the exterior of such black hole spacetimes. However, it does not constrain the parameter lo anyway. 

Thus, it was argued in e.g. [41] that the entropy of such black holes should be ^irpo A, where A 
is the horizon area and po is constant (corresponding to some finite value of lo), i.e not necessarily 
it equals p^°'> = 

However, when we consider the de Sitter black hole spacetimes, it turns out that p is not only 
a constant between the black hole and the cosmological event horizon, but also we must have 
uj = oo |42j . In other words, for such spacetimes, we must have po = p^^^ = It is then clear 

that the total entropy of the de Sitter black holes in Brans-Dicke theory equals m- Clearly, this 
is qualitatively different from the asymptotically flat black hole spacetimes. 

5 Discussions 

In this work we have considered the thermodynamics of stationary axisymmetric de Sitter black hole 
spacetimes. We have utilized the formalism developed using the Gibbons-Hawking-York surface 
couterterm and near horizon symmetries [22]-|25j to derive the total entropy of such two horizon 
spacetimes. To the best of our knowledge, this has not been done before. Since the spacetimes we 
have considered are endowed with two Killing horizons, the surface counterterm method provides 
us a natural convenience to deal with them. We have used a very general geometric framework to 
perform the derivation of the entropy in Sec.s 2 and 3. We have also seen in Sec. 3 that in order to 
do that, we have to choose uniquely, constants in the mode functions preserving the near horizons’ 
geometries. Such restriction is not present for single horizon spacetimes. 

This entire analysis can be thought of as doubly-local, instead of local as the single horizon 
spacetimes. Since nowhere in our analysis we assumed any precise asymptotic behaviour of the 
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metric, our analysis is absolutely valid for the Nariai class de Sitter black holes, where the black 
hole and the cosmological event horizons have comparable sizes. 

We also emphasize here that considering two boundary integrals in such two horizon spacetimes 
is the most natural choice. 

After performing this general analysis, we have considered known non-trivial stationary axisym- 
metric solutions belonging to the Kerr-Newmann- or more general Plebanski-Demianski-de Sitter 
classes, and have demonstrated that for all of them the general derivation of the total entropy holds. 
We have also considered two cases of non-minimal couplings. In particular, we have pointed out 
a qualitative difference of the de Sitter black hole entropy in the presence of a Brans-Dicke scalar 
field, compared to the asymptotically flat spacetimes. 

The analysis we have done above holds well as long as the spacetime is not exactly extremal, 
as the very use of Eq. (fT^ requires, even though k could be ‘small’, but could never be vanishing. 
This is just because of the fact that the extremal black holes are qualitatively different objects from 
the usual non-extremal or even near-extremal ones. The vanishing of the surface gravity indicates 
vanishing temperature and this leads to a debate about the entropy of such black holes : should the 
entropy of such black holes vanish? We refer our reader to |44] for a recent discussion and a novel 
proposal about the calculation of entropy of extremal Schwarzschild-de Sitter spacetime, where the 
black hole and cosmological event horizon coincide. In that spacetime, there exists no timelike 
Killing vector field. It was shown in [53] that there exists a vector field, which is not Killing on the 
bulk, but becomes so and also null on the coinciding horizons, thereby enhancing the symmetry of 
the spacetime near the horizons. This is in accordance with the usual requirement that the near 
horizon structure must contain the conformal group. Rather curiously, this particular vector field 
gives bifurcate Killing horizons, spatially coincident with the original horizons we started with, 
which might give rise to a nonvanishing entropy. It should be interesting to attempt to relate the 
current formalism we have used, with that of [44] . Precisely, for non-extremal spacetimes, where 
the black hole and the cosmological horizon are separated, one might look for such a vector field 
and check its near horizon properties. The enhanced algebra containing this vector field might yield 
meaningful result for all values of k, and hence might give a meaningful notion of the yet not-well 
understood k —>■ 0 limit. This procedure seems to have some qualitative similarity with what we 
did for the stationary axisymmetric spacetimes - we found a specific timelike non-Killing vector 
held x“ in the bulk, becoming Killing and null on both the horizons. 

As we also have emphasized in the beginning, the current work is a step towards understanding 
aspects of de Sitter black hole thermodynamics when we treat the entire two horizon spacetime as 
a thermodynamic system as a whole. The variation of the total entropy of such spacetimes gives a 
Smarr formula predicting an effective equilibrium temperature. Thus the next step is to understand 
how or under what circumstances or with what choices of the field mode functions one can actually 
derive such thermal states using field theory. We hope to return to this issue in future works. 


Acknowledgement 

This research was implemented under the “ARISTEIA 11” Action of the Operational Program 
“Education and Lifelong Learning” and is co-funded by the European Social Fund (ESF) and 
Greek National Resources, when I was a post doctoral researcher at ITCP and Dept, of Physics of 
University of Crete, Heraklion, Greece. My current research is supported by lUCAA, Pune, India. I 
thank Amit Ghosh, Avirup Ghosh and Amitabha Lahiri for useful discussions and encouargement. 


References 

[1] G. W. Gibbons and S. W. Hawking, Phys. Rev. D 15, 2738 (1977). 

[2] D. Kastor and J. H. Traschen, Phys. Rev. D 47, 5370 (1993). 


15 


[3] N. Goheer, M. Kleban and L. Susskind, JHEP 0307, 056 (2003). 

[4] M. I. Park, Class. Quant. Grav. 26, 075023 (2009). 

[5] H. Saida, Prog. Tlieor. Pliys. 122, 1515 (2010). 

[6] H. Saida, Prog. Theor. Phys. 122, 1239 (2010). 

[7] M. Urano, A. Tomimatsu and H. Saida, Class. Quant. Grav. 26, 105010 (2009). 

[8] D. Kastor, S. Ray, J. Traschen, Class. Quantum Gravity 26, 195011 (2009). 

[9] B. P. Dolan, Class. Quantum Gravity 28, 125020 (2011). 

[10] B. P. Dolan, D. Kastor, D. Kubiznak, R. B. Mann, J. Traschen, Phys. Rev. D 87, 104017 
(2013). 

[11] S. Bhattacharya and A. Lahiri, Eur. Phys. J. C 73, 2673 (2013). 

[12] S. Carlip, Class. Quant. Grav. 15, 3609 (1998). 

[13] S. Das, A. Ghosh and P. Mitra, Phys. Rev. D 63, 024023 (2001). 

[14] O. Dreyer, A. Ghosh and J. Wisniewski, Class. Quant. Grav. 18, 1929 (2001). 

[15] S. Silva, Class. Quant. Grav. 19, 3947 (2002). 

[16] S. Carlip, Gen. Rel. Grav. 39, 1519 (2007). [Int. J. Mod. Phys. D 17, 659 (2008)] 

[17] O. Dreyer, A. Ghosh and A. Ghosh, Phys. Rev. D 89, 024035 (2014). 

[18] G. Kang, J. i. Koga and M. I. Park, Phys. Rev. D 70, 024005 (2004). 

[19] J. L. Cardy, Nucl. Phys. B 270, 186 (1986). 

[20] H. W. J. Bloete, J. L. Cardy and M. P. Nightingale, Phys. Rev. Lett. 56, 742 (1986). 

[21] S. Carlip, Int. J. Mod. Phys. D 23, 1430023 (2014). 

[22] B. R. Majhi and T. Padmanabhan, Phys. Rev. D 85, 084040 (2012). 

[23] B. R. Majhi and T. Padmanabhan, Phys. Rev. D 86, 101501 (2012). 

[24] B. R. Majhi and T. Padmanabhan, Eur. Phys. J. C 73, 2651 (2013). 

[25] B. R. Majhi, Phys. Rev. D 90, no. 4, 044020 (2014). 

[26] B. R. Majhi, JCAP 1405, 014 (2014). 

[27] K. Meng, Z. N. Hu and L. Zhao, Class. Quant. Grav. 31, 085017 (2014). 

[28] P. C. W. Davies and T. M. Davis, Found. Phys. 32, 1877 (2002). 

[29] K. Maeda, T. Koike, M. Narita and A. Ishibashi, Phys. Rev. D 57, 3503 (1998). 

[30] S. Shankaranarayanan, Phys. Rev. D 67, 084026 (2003). 

[31] F. L. Lin and Y. S. Wu, Phys. Lett. B 453, 222 (1999). 

[32] B. Wang, E. Abdalla and R. K. Su, Phys. Lett. B 503, 394 (2001). 

[33] L. C. Zhang, M. S. Ma, H. H. Zhao and R. Zhao, Eur. Phys. J. C 74, 3052 (2014). 

[34] S. Bhattacharya, Phys. Rev. D 88, 044053 (2013). 

[35] R. M. Wald, “General Relativity,” Chicago Univ. Pr. (1984)- 

[36] Z. W. S. Chong, M. Cvetic, H. Lu and C. N. Pope, Phys. Rev. Lett. 95, 161301 (2005). 

[37] J. B. Griffiths and J. Podolsky, Int. J. Mod. Phys. D 15, 335 (2006). 

[38] G. W. Gibbons, H. Lu, D. N. Page, and C. N. Pope, Phys. Rev. Lett. 93, 171102 (2004). 

[39] G. W. Gibbons, H. Lu, D. N. Page, and C. N. Pope, J. Geom. Phys. 53, 49 (2005). 

[40] S. Bhattacharya and H. Maeda, Phys. Rev. D 89, 087501 (2014). 

[41] V. Faroni, Entropy 12, 1246 (2010). 

[42] S. Bhattacharya, K. F. Dialektopoulos, A. E. Romano and T. N. Tomaras, Phys. Rev. Lett. 
115, no. 18, 181104 (2015). 

[43] S. W. Hawking, Commun. Math. Phys. 25, 167 (1972). 

[44] S. Stotyn, Can. J. Phys. 93, no. 9, 995 (2015), [arXiv:1502.02737 [gr-qc]]. 


16 


